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ABSTRACT 

We consider the radial structure of radiatively inefficient hydrodynamic accretion 
flows around black holes. We show that low-viscosity flows consist of two zones: an 
outer convection-dominated zone and an inner advection-dominated zone. The transition 
between these two zones occurs at ~ 50 Schwarzschild radii. 
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1. Introduction 

Observations of low luminosity black holes in galac- 
tic nuclei and X-ray binaries have been modeled 
within the framework of a radiatively inefficient accre- 
tion flow (RIAF). In these models, radiative losses are 
small because of the low particle density of the accret- 
ing plasma at low mass accretion rates (see Narayan, 
Mahadevan & Quataert 1998 and Kato, Fukue & Mi- 
neshige 1998 for reviews). 

Hydrodynamical simulations of low-viscosity RI- 
AFs in two dimensions (Igumenshchev & Abramow- 
icz 1999, 2000; Stone, Pringle & Begelman 1999) 
and three dimensions (Igumenshchev, Abramowicz & 
Narayan 2000) have shown that these flows are con- 
vectively unstable (sec also Narayan & Yi 1994, 1995), 
and that the convection has a significant effect on 
the structure of the flow. Narayan, Igumenshchev 
& Abramowicz (2000, hereafter NIA) and Quataert 
& Gruzinov (2000, hereafter QG) constructed self- 
similar solutions for convection-dominated accretion 
flows (CDAFs), and used these to explore the relevant 
physics and to explain the results of the numerical 
simulations. 

In this paper we show via analytical arguments 
that the radial structure of a low-viscosity RIAF 
consists of two distinct zones: an outer convection- 
dominated zone and an inner advection-dominated 
zone (§2). We present global solutions of an idealized 
height-integrated set of equations (§3) and show that 
the radial structure of the advection-dominated zone 
differs significantly from the canonical self-similar 
solution for an advection-dominated accretion flow 
(ADAF, Narayan & Yi 1994, 1995). This is because 
of the strong influence of the boundary conditions. 
We then compare the theoretical results with numeri- 
cal simulations (§4) and find that there is good agree- 
ment. In particular, both the global solutions and the 
numerical simulations show that the transition from 
the inner ADAF zone to the outer CDAF zone occurs 
at a radius Rac ~ bQRg, where Rg = 2GM/c^ is 
the gravitational radius of the accreting black hole of 
mass M. 



2. Analysis of Self-Similar CDAFs 

We consider a rotating accretion flow with a self- 
similar scaling of the height-averaged angular velocity 
Q. and isothermal sound speed Cg, and the scale height 
H: 



ctiR) = 4vj, <x R-\ (1) 

H{R) = Cs/Qk = coR, 

where VIk — \/GMJR? = vk/R is the Keplerian 
angular frequency at radius R, and Q.q and cq are di- 
mensionless constants. We consider a power-law ra- 
dial dependence for the height averaged density, 



p{R) oc R-°-, 



(2) 



which corresponds to a radial velocity varying as 

VR (R) = M /AttRHp oc (3) 

where M is the mass accretion rate. The index a is 
equal to 1/2 for a self-similar CDAF (NIA) and 3/2 
for a self-similar ADAF (Gilham 1981; Begelman & 
Meier 1982; Narayan & Yi 1994, 1995). We refer the 
reader to NIA for the calculation of the constants Oq 
and Cq. 

Here we focus on the consistency of the self-similar 
CDAF solution. For this purpose we consider the 
internal energy equation describing a RIAF, 



P^'^^dR + R^dR^'^^^^ 



radi 



(4) 



where s is the specific entropy, is the outwardly 
directed convective energy flux, Q'^ is the rate of en- 
ergy dissipation per unit volume, and Qrad is the rate 
of energy loss per unit volume through radiation. By 
deflnition, radiative losses are small in a RIAF, so we 
neglect Qrad in equation (4). The first term on the left 
hand side of equation (4) describes the inward advec- 
tion of the internal energy. For the scaling laws given 
in (1) and (2) this advection term takes the form. 



ds_ ^ Mc^ clRg 
^"""^ dR AttRH 2 i?2 



1 



7-1 



— a 



(5) 



where c is the speed of light and 7 is the adiabatic 
index of the accreting gas. For a rotating accretion 
flow, the energy dissipation rate Q"*" can be written 
as 

Q+ = -{Jv + Jc)^, (6) 

where Jy and Jc are the angular momentum fluxes due 
to viscosity and convective motions, respectively. In 
the case of an ADAF (a = 3/2), convection is ignored 
and we set Jc = = 0. Since viscosity transports 
angular momentum outwards, i.e. J„ > 0, the dissi- 
pation rate is a positive quantity, and in equation 
(4), (3+ balances the advection term (5). 
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The situation is very different in a CDAF. NIA 
considered a self-similar CDAF in the limit of a non- 
accreting "convective envelope" , in which vji ^ and 
the advection term (5) vanishes. How can a viscous 
differentially-rotating fluid configuration not accrete? 
This is possible because, as demonstrated explicitly 
via three-dimensional numerical simulations (Igumen- 
shchev et al. 2000) and earlier through analytic meth- 
ods (Ryu & Goodman 1992; Stone & Balbus 1996), 
axisymmetric convection transports angular momen- 
tum inwards, i.e. Jc < 0, whereas viscosity transports 
it outwards. Thus, by ensuring that convection and 
viscosity cancel each other exactly, i.e. Jy + Jc — 0, 
the convective envelope is able to achieve a configu- 
ration with no net angular momentum transport and 
hence no mass accretion. Since Jy + Jc ^ 0, an ad- 
ditional consequence is that = (see eq. [6]). 
The energy equation (4) then has only one term left, 
namely the divergence of the convective flux, and the 
only way to satisfy this equation is for the divergence 
of the flux to vanish, i.e. to have Fc{R) oc R~^. 

The non-zero outward energy flux Fc corresponds 
to a "convective luminosity" , Lc = 4ttRHFc, which is 
independent of radius. The source of this luminosity 
is not specified in the self-similar solution and is for- 
mally located at R — 0. NIA argued that in the case 
of a real CDAF there would be a small but non-zero 
accretion rate, M ^ 0, and the convective luminos- 
ity would be supplied by the gravitational release of 
binding energy in the innermost region of the accre- 
tion flow. Thus, we may write 

ic-eA/c^ (7) 

where e is the "convective efficiency" . A non-zero ac- 
cretion rate requires a non-zero net angular momen- 
tum flux directed outward, Jy + Jc > 0, so according 
to (6) one would also have > 0.| 

We now estimate the importance of different terms 
in the energy equation (4) . For this we integrate equa- 
tion (4) over a spherical volume from R to infinity. 
After the integration, the terms which correspond to 
the advection term (5) and the dissipation term (6) 
both scale oc R~^, whereas the term describing the 
convective energy transport is equal to Lc — con- 
stant. Thus, convective energy transport dominates 
the physics at large radii, and the self-similar CDAF 

^Taking into account that the release of gravitational and rota- 
tional energy is the ultimate source of any dissipation in the 
accretion flow, one can conclude oc {GM/R)pvfi/R oc _R~*. 



solution is a good approximation for describing this 
region of a RIAF. At small radii, however, the advec- 
tion term and the viscous dissipation term become 
more important. By comparing the inward advective 
energy flux, 

^ /; i,RH .„„r|l .fl 4 (-1^ - a) §M.^ 

(8) 

with the convective luminosity Lc (eq. [7]), we may 
determine the transition radius Rac that separates 
the outer convection-dominated zone and the inner 
advection-dominated zone. Using the estimate of Cq 
from QG or NIA's Appendix, and applying the con- 
dition Eadv ~ Lc, we obtain 

Rac 1 3-7 

Rg ^ e2{j + 3)- ^ ^ 

For a quantitative estimate of Rac we need to deter- 
mine the value of the convective efficiency parameter 
s. Since it is hard to determine this from first prin- 
ciples, we must either solve for e as an eigenvalue of 
the global problem (§3) or determine it from numeri- 
cal simulations (§4). 

3. Global CDAFs 

In this section we present a sample height inte- 
grated global CDAF model which satisfies proper 
inner and outer boundary conditions. This is the 
CDAF analogue of the global ADAF models calcu- 
lated by Narayan, Kato, & Honma (1997) and Chen, 
Abramowicz, & Lasota (1997). The solution sheds 
further light on the CDAF- ADAF transition discussed 
in §2. 

In the absence of radiation, the conservation of 
mass, angular momentum, and total energy can be 
expressed via conserved fluxes: 

M = -AttRHpvr = const, (10) 

At:Fl = -MR^VL + AirRH^Jy + Jc) = const, (11) 

47rFE = -MB+47rRHn{Jy + Jc)+ATrRHFc = const, 

(12) 

where B = 0.5w^ + (f> + O.bR^fl^ + 70^/(7 - 1) is 
the Bernoulli function, and we take the Paczyhski & 
Wiita (1980) form of the gravitational potential, (j) = 
—GM I {R — Rg), to mimic a Schwarzschild black hole. 
Note that the radial derivative of equation (O) yields 
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equation (Q) with Qrad = 0. In addition to equations 
(10)-(E2|), we need the radial momentum equation: 



dvB. „„2 1 dp 

''^llR='^^ -dR--pdR^ 

where p = pel is the pressure. 

We use the expressions for Jy, Jc and Fc from 
NIA's Appendix except that we employ a causal pre- 
scription for the convective transport of energy and 
angular momentum.|^ We take the convective diffu- 
sion coefficient to be oc (1 — v^/cD"^ for < Cc and 
otherwise (cf. Popham & Narayan 1992; Narayan 
1992), where we choose Cc — aCg (i.e. roughly the 
convective velocity). The condition \vr\ — Cc defines 
a transition radius, RaCi inside of which there is no 
convective transport of angular momentum and en- 
ergy. For R < Rac , the inflow velocity of the gas ex- 
ceeds the characteristic convective velocity, and con- 
vective transport is less important than bulk advec- 
tive transport (the convective transport vanishes in 
our simplified model). Outside of Rac the inflow ve- 
locity of the gas is less than the characteristic con- 
vective velocity and the flow is roughly "static" , as in 
the self-similar CDAF solution. 

We note that Jy depends on the value of the 
Shakura-Sunyaev viscosity parameter a, so the global 
solution depends on the choice we make for this pa- 
rameter. Other parameters are the adiabatic index 7 
and the outer radius of the accretion flow Rout ■ In ad- 
dition, although the causal prescription for convective 
transport is well motivated (e.g., Narayan 1992), the 
proper cutoff velocity (i.e., Cc) is somewhat uncertain. 

Equations ([l^l-dll) 

are four first order differen- 
tial equations. Since we are interested in solutions in 
which the gas goes through a sonic point on its way 
into the black hole, the sonic radius Rs is an addi- 
tional unknown which needs to be determined. The 
3 conserved fluxes, M, Fi^, and Fe are also eigen- 
values determined in a self-consistent manner. This 
implies that the convective efficiency e of equation (0) 
is an eigenvalue of the global model, rather than a free 
parameter as in the self-similar solution. The transi- 
tion radius Rac is a final eigenvalue. For R < Rac, 
Fc = and Jc = because of the causal prescrip- 
tion for convective transport. One can show that the 



vanishing of Fc means that Fe and Fl are not inde- 
pendent, Fl — Fe{R^^/ B)\R=ii^. Thus, we require a 
total of 8 boundary conditions to solve for the 4 vari- 
ables (p, Ufl,, Cs, ri) and 4 independent eigenvalues 
{Rs, Rac, M, Fe). 

At the outer radius of the fiow we specify p, Cs, 
and fl. The choice of p is arbitrary and simply sets 
the scaling for the density in the problem. We use 
the self-similar CDAF values for Cs and from QG 
at the outer boundary. Igumenshchev & Abramow- 
icz (1999, 2000) showed from numerical simulations of 
RIAFs that strong convection is present only for rela- 
tively small values of the viscosity parameter, a < 0.1. 
For such values of a, Rac > Rs- At the sonic point 
R = Rs, convection is thus unimportant and the 2 
boundary conditions at this radius are identical to 
those in the global ADAF models of Narayan et al. 
(1997). We also apply the no-torque condition at the 
sonic point The final two boundary conditions are 
at the transition radius Rac- We guess a value for 
Rac and find the solution from the outer boundary to 
Rac and from Rac to the inner boundary (the sonic 
point), requiring = Cc at Rac- We then iterate 
until all flow variables are continuous across Rac- 

Figure 1 shows the density, radial velocity, and con- 
vective velocity Cc for a global solution with a — 0.03 
and 7 = 1.5. The sonic point of the solution is at 
R w 2.6Rg, while the transition radius is located at 
Rac ~ 35i?g. The angular momentum flux in the so- 
lution is AttFl ~ —l.GAMRgC while the energy flux is 
4:1tFe ~ 0.0045Mc'^. At large radii the energy flux is 
entirely carried by convection, so this corresponds to 
a convective efficiency of e « 0.0045. For this value 
of e and 7 = 1.5, equation (9) predicts Rac = 'iTRg, 
which is very close to the numerically determined lo- 
cation of the transition radius. 

The density profile in Figure 1 is reasonably well 
described by the p oc R^-^^^ power law of the self- 
similar CDAF solution at all radii. For R > Rac con- 
vective transport indeed dominates the energy trans- 
port and the scaling p oc R^^^^ is needed in order 
to carry a constant convective luminosity outwards 
(see §2 and NIA, QG); in this regime the self-similar 
CDAF solution describes the flow structure quite well. 

For R < Rac, however, convective transport of en- 
ergy and angular momentum vanishes in our causal 



^For numerical reasons, we were unable to find global solutions 
without employing a causal transport prescription. This is 
somewhat unsatisfying because the causal prescription explic- 
itly imposes an inner ADAF zone on the global CDAF solution. 



^Narayan et al. (1997) applied the no-torque condition close to 
the black hole horizon; we do not expect this to make a large 
difference. 
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prescription. Hero, the radial inflow velocity exceeds 
the convective velocity and bulk (advective) transport 
of energy dominates. The density profile is nonethe- 
less close to the CDAF scaling of R~^l'^ and much 
flatter than the Br^l"^ profile of a self-similar ADAF. 
The reason for this is that the gas is rapidly acceler- 
ating toward the sonic point on its way into the black 
hole. Quite generally, this rapid acceleration leads to 
a steep radial velocity profile and thus a flat density 
profile (see, e.g., Narayan ct al.'s 1997 and Chen et 
al.'s 1997 global ADAF models for similar results). 

4. Comparison with Numerical Simulations 

Igumenshchev (2000) and Igumenshchev & Abramow- 
icz (2000) found that the convective efficiency is e ~ 
0.003 — 0.01 in two-dimensional simulations of RI- 
AFs, with only a weak dependence on the viscosity 
parameter a and the adiabatic index 7. We expect 
that three-dimensional models of convective RIAFs 
will show values of e in the same range, because of 
the close qualitative and quantitative similarities be- 
tween the numerical results in two and three dimen- 
sions (Igumenshchev et al. 2000). Taking a charac- 
teristic value of e = 0.005, we obtain from equation 
(9): Rac ^ 30i?,g for 7 = 5/3 and Rac ^ 40i?g for 
7 = 4/3. These values for e and Rac are in good 
agreement with the values obtained with the global 
model described in §3. We now try to estimate Rac 
directly from the numerical simulations. 

The presence of two distinct zones in a RIAF, 
namely an outer convection-dominated zone and an 
inner advection-dominated zone, may be demonstrated 
in a numerical simulation by considering the "mass 
inflow rate" Mi„ as a function of radius. We de- 
fine Min{R) = J pvRdS, where the surface integral 
is limited to those gas elements at radius R that have 
vr < 0. Because a CDAF is convectively turbulent 
and accretes slowly, roughly half the mass at any ra- 
dius at any time will be flowing in and half will be 
flowing out. This gives Mi„ ~ AttRH pv/2, where v 
is the rms velocity of the turbulent eddies. It c:an be 
shown that v ^ av^ oc (NIA). Therefore, for 

p oc and H (x R, one finds that Mm oc R. 

In an ADAF on the other hand, convection is absent 
and so the mass inflow rate is independent of radius, 
Mi„{R) = const. 

Figure 2 shows Min{R) (solid line) for a two- 
dimensional simulation of a RIAF with a = 0.01 
and 7 = 5/3 (the numerical technique is similar to 



that used by Igumenshchev & Abramowicz 2000). 
This simulation was done with the Paczynski & Wi- 
ita (1980) pseudo-potential. The plot clearly shows 
that the accretion flow consists of two different zones: 
an inner advection-dominated zone, in which Mj„ = 
const, and an outer convection-dominated zone, in 
which Min oc R. The transition radius between the 
two zones is located at Rac — 50 i?g, which is roughly 
consistent with the scaling given in equation (9). 

Figure 3 shows a snapshot of the velocity stream- 
lines of the particular numerical simulation analyzed 
in Fig. 2. We see numerous vortices associated with 
the convective motions in the flow. Interestingly, un- 
like in the oversimplified global model described in 
§3, we see that convective motions are present even in 
the inner regions of the RIAF {R. < Rac), despite the 
dominance of advection here. In fact, the convective 
blobs primarily originate in the advection-dominated 
part and move outward. According to our estimate 
(see equation [5] and [7]), the major energy release 
happens in the advection-dominated part of RIAFs, 
where the advection and dissipation terms are almost 
equally important and dominate the convective en- 
ergy transport. Only a small relative excess of the 
energy dissipation rate over the inward advection of 
energy supports turbulence in RIAFs and provides 
the convective luminosity Lc- This excess determines 
the value of the convective efficiency e. 

5. Discussion and Conclusions 

We have considered the radial structure of radia- 
tively inefficient accretion fiows (RIAFs). Based on 
an analysis of self-similar solutions we conclude that 
the flow consist of two parts: an outer convection- 
dominated part where the radial velocity is highly 
subsonic and an inner advection-dominated part where 
the gas flows rapidly inwards. The location of the 
transition radius Rac between these two parts de- 
pends primarily on the convective efficiency e (see 
equation [9]). We estimate e « 0.003 — 0.01 and 
Rac ~ 50Rg directly from numerical simulations 
(§4). In addition, in ID steady state global cal- 
culations of the structure of RIAFs, which include 
proper boundary conditions, e and Rac can be de- 
termined as eigenvalues. The global calculations yield 
Rac ~ 35i?g and e w 0.0045 (§3), in good agreement 
with the simulations. 

Convective motions and the associated outward en- 
ergy flux originate in the inner advection-dominated 
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part of the flow. This energy flux supports the tur- 
bulent structure of the flow in the outer convection- 
dominated part. The radial profiles of density and 
velocity in the inner advection-dominated part of the 
flow differ significantly from self-similar scalings due 
to the influence of boundary conditions (§3). 

Our present results are based on a viscous hydro- 
dynamical (HD) approach to studying RIAFs. The 
HD approach is not self-consistent because it requires 
one to assume the presence of viscosity with an a pri- 
ori unknown strength. A fully consistent approach 
can be provided within the framework of magneto- 
hydrodynamics (MHD), in which torques due to mag- 
netic interactions transport angular momentum. It 
is not presently clear how the concept of CDAFs ob- 
tainc;(l in the HD approach will be modified after con- 
struction of MHD models of RIAFs. 

We do not expect significant qualitative differences 
between HD and MHD models if the gravitational and 
rotational energies are released locally in the accretion 
flow, e.g., by magnetic reconnection and/or a turbu- 
lent energy cascade. If, however, the free energy is 
removed from the bulk of the accretion flow via non- 
local interactions, such as the stretching of large-scale 
magnetic field lines resulting in an efficient Poynting 
flux, or the formation of MHD driven winds/outflows, 
one would expect qualitative differences between HD 
and MHD models. In addition, one would expect a 
significant modification of the CDAF solution if an- 
gular momentum is efficiently removed from the ac- 
cretion flow by a large-scale magnetic field in a disk 
corona (contrary to the CDAF model in which angu- 
lar momentum transport is due to local viscous inter- 
actions). 

MHD simulations of spherical accretion flows (Igu- 
menshchev &: Narayan 2001) support the idea of lo- 
cal generation of energy via magnetic reconnection, 
which leads to the development of efficient convec- 
tion. The structure of these models closely resem- 
bles CDAFs. Recent MHD simulations of rotating 
flows do show unstable motions, which leads to the 
development of efficient turbulence, but they do not 
show clear behaviour to either support or contra- 
dict the CDAF model. For example, Machida. Mat- 
sumoto, & Mineshige (2001) claim that the structure 
of their MHD models closely resembles that found 
in HD CDAFs. However, Hawley, Balbus, & Stone 
(2001) stress that their models do not demonstrate 
convective motions and that the models are instead 
turbulent due to the magneto-rotational instability in 



shearing flows. This discrepancy indicates that fur- 
ther development of the MHD models is required. 

Recently, Balbus (2001) has discussed the stabil- 
ity criteria for rotating magnetized accretion flows. 
These criteria replace the H0iland criteria for a plasma 
with thermal conduction along the magnetic field 
lines. According to Balbus's criteria a magnetized 
fluid is unstable even if it is marginally stable to the 
H0iland criteria. The new criteria have been obtained 
in the short-wavelength approximation, and it is not 
straightforward to know how the instability works on 
larger scales, which might be more important for the 
structure of an accretion flow. It is also interesting 
to clarify what Balbus's criteria imply for axisym- 
metric convection in the presence of a predominantly 
toroidal magnetic field. Such convection is expected 
to be present in magnetized CDAFs. 

Assuming that the radial structure of RIAFs de- 
scribed in the present paper is basically correct, it is 
interesting to ask what the implications are for ob- 
servations of accreting black holes. It has been pro- 
posed that for luminosities below a few percent of the 
Eddington luminosity, accretion flows around black 
holes switch to a RIAF (cf. Narayan et al. 1998). 
In the scenario proposed by Narayan (1996) and de- 
veloped by Esin, McChntock & Narayan (1997), the 
accretion flow consists of a radiatively efficient thin 
disk outside a transition radius Rtr and an ADAF in- 
side Rtr- Spectra computed with such a model agree 
well with observations of black hole X-ray binaries in 
the "intermediate state" and "low state" (Esin et al. 
1997, 1998). The models typically require Rtr of or- 
der lORg to a few tens of Rg. In a recent study, Esin 
et al. (2001) determined that the black hole binary 
XTE J1118+480 had Rtr 50i?g at a time when its 
X-ray luminosity was ^ 10^^ LEdd- According to the 
present paper, an ADAF is a good approximation to 
the structure of a RIAF so long as i? < 50i?g. There- 
fore, the above models are consistent. 

For yet lower luminosities, e.g., black holes in the 
"quiescent state" , the transition radius in models is 
typically at Rtr ~ 10^ - 10"^ Rg, and in some models 
there is no thin disk at all (see Narayan et al. 1998 for 
examples). Published models of the qiuescent state 
have generally made use of a pure ADAF (or a spher- 
ical Bondi flow, e.g. Melia 1992), whereas our present 
results suggest that the flow might take the form of 
a CDAF for R ^ bORg. How will this affect the pre- 
dicted spectrum? The principal difference between 
an ADAF (or a Bondi flow) and a CDAF is the value 
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of the index a (cf cq. [2]) in the radial density pro- 
file: a = 3/2 for a self-similar ADAF, a = 1/2 for a 
self-similar CDAF. Quataert & Narayan (1999) stud- 
ied models with different values of a and showed that 
these models can produce similar spectra provided the 
fraction of viscous heat energy that goes into electrons 
(which is assumed to be close to zero for an ADAF) 
is appropriately adjusted. For a CDAF, perhaps as 
much as half the heat energy will need to go into the 
electrons to match observations (e.g. Ball, Narayan 
& Quataert 2001). Detailed models are awaited. 
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Fig. 1.— Global CDAF model with a = 0.03 and 
7 = 1.5. The density profile p is in arbitrary units. 
Inside R = Rac ~ 35Rg the inflow speed exceeds the 
local convective velocity (Cc) and energy transport in 
the flow is dominated by inward advection. At larger 
radii, the flow is roughly static {\vr\ ^ Cc) and energy 
is transported outward by convection. 
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Fig. 2. — Instantaneous mass inflow rate Mm (solid 
line) as a function of radius i? in a 2D hydrodynam- 
ical model of a radiatively inefficient accretion flow 
with 7 = 5/3 and a = 0.01. The values of Mi„ 
are calculated by adding all the inflowing gas ele- 
ments at a given R. Min is normalized to the net 
accretion rate Mq. Except near the outer boundary, 
R ^ 2 X 10^ Rg and in the inner region, R <^ 10^ Rg, 

the profile Min{R) shows good agreement with the 
prediction of the self-similar solution for convection- 
dominated accretion flows, Mi„ oc i?, shown by the 
long-dashed line. Deviation of Mi„ from the self- 
similar CDAF scaling in the inner region is due to 
the strong inward advection of internal energy, as dis- 
cussed in the text. 
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Fig. 3. — Snapshot of streamlines from a 2D hydro- 
dynamical model of a radiatively inefficient accretion 
flow with a = 0.01 and 7 = 5/3. A projection of 
streamlines on the meridional cross-section is shown, 
with the black hole located at the origin. The flow 
pattern is highly time dependent and consists of nu- 
merous temporal vortices of different spatial scales. In 
the innermost region of the flow the energy balance 
is dominated by inward advection of internal energy. 
Convectivc blobs originate in this region and move to- 
ward the outer convection-dominated region at radii 
^ lO^iJg. 
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